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1.  Introduction 


Guided  Lamb  wave  nondestructive  evaluation  is  a  tool  used  in  structural  health  monitoring  of 
aircraft,  civil,  and  mechanical  infrastructures.  The  goal  is  to  detennine  how  wave  velocity  and 
mode  shape  are  influenced  by  changes  in  geometry  and/or  boundary  conditions  caused  by 
structural  damage  or  degradation.  These  waves  propagate  in  thin  plates  and  plate-like  structures 
and  are  formed  by  the  interference  of  multiple  reflections  and  mode  conversions  of  longitudinal 
and  shear  waves  at  the  plate  surfaces.  Guided  waves  are  typically  generated  and  detected  using 
conventional  piezoelectric  transducers,  oriented  either  flat  or  at  an  angle  with  respect  to  the 
surface.  Two  types  of  waves  propagate  a  symmetric  wave  and  an  antisymmetric  wave,  each  with 
multiple  modes  and  at  different  frequency  and  speeds  that  complicate  data  analysis.  To  avoid 
analysis  complications,  fundamental  symmetric  (So)  and  antisymmetric  (A0)  modes  at  low- 
frequency  thickness  values  are  typically  used  for  inspections  by  limiting  source  transducer 
bandwidth.  Processing  and  analysis  of  guided  wave  signals  are  typically  performed  in  the  time 
and  frequency  domain.  Other  popular  analysis  methods  include  integrated  time-frequency 
domain  and  wavelet  analysis.  In  some  cases,  baseline  signals  are  subtracted  from  inspection 
signals  to  enhance  results. 

A  wide  range  of  work  has  been  reported  on  the  interaction  of  guided  waves  (Lamb  waves)  with 
damage.  In  these  works,  the  majority  of  analysis  was  done  through  examining  transmission  and 
reflection  coefficients.  For  example,  Lamb  wave  interaction  with  straight  sided  notches  was 
investigated  by  Alleyne  and  Cawley  (1 )  through  a  finite-element  model  and  two-dimensional 
(2-D)  Fourier  analysis  of  low-order  modes.  Their  results  showed  that  the  reflection  and 
transmission  coefficients  of  individual  Lamb  wave  modes  to  notch  geometry  and  orientation  are 
dependent  on  frequency,  plate,  and  wave  mode.  Lowe  et  al.  (2)  and  Lowe  and  Diligent  (3)  also 
investigated  Lamb  wave  interaction  with  a  rectangular  notch  for  first-order  symmetric  and 
antisymmetric  modes.  They  used  reflection  coefficients  to  quantify  results  that  showed  a  sharp 
rise  in  reflection  coefficient  for  a  notch  depth-to-plate  thickness  of  50%  and  greater. 

Benmeddour  et  al.  ( 4 )  used  a  method  to  predict  the  propagation  of  fundamental  Lamb  wave 
modes  in  an  isotropic  structure  containing  symmetric  notches;  power  reflection  and  transmission 
coefficients  were  used  in  his  analysis.  Differential  guided  wave  features  were  used  by  Michaels 
et  al.  (5)  to  classify  structural  changes  such  as  holes  and  temperature  variations.  Classifiers 
utilizing  time  and  frequency  domain  features  were  compared  to  classifications  based  on  time- 
frequency  representations.  Both  classification  schemes  were  able  to  discriminate  between 
environmental  and  structural  changes  in  a  small  aluminum  plate.  In  this  report,  guided  waves 
transmitted  across  a  notch  (simulated  damage)  in  a  metallic  plate  are  numerically  generated  and 
examined  with  various  analysis  methods  to  detennine  the  wave  defect  sensitivity  and  potential 
damage  characterization. 
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2.  Fundamental  of  Lamb  Waves 


Guided  waves  or  Lamb  waves  are  formed  by  the  interference  of  multiple  reflections  and  mode 
conversions  of  longitudinal  and  shear  waves  at  the  free  surfaces  of  a  plate.  These  waves  are 
typically  generated  and  detected  using  conventional  piezoelectric  transducers,  oriented  either  flat 
or  at  an  angle  with  respect  to  the  surface,  and  are  used  to  detect  defects  and  measure  elastic 
properties  of  thin  isotropic  materials  and  laminated  composite  plates.  Two  types  of  waves 
propagate  a  symmetric  wave  and  an  antisymmetric  wave.  Each  of  these  wave  types  propagates 
with  multiple  modes. 

By  solving  a  boundary  value  problem  for  a  free  plate  or  by  considering  the  reflection  of  waves  at 
plate  boundaries,  the  Rayleigh-Lamb  frequency  equations  (dispersion  equations)  can  be  found 
(6-9).  For  a  unifonn  traction-free  isotropic  plate,  the  equation  for  symmetric  modes  is 

tan  (<?/?)  _ 
tan  ( ph j 

For  antisymmetric  modes, 

tan (qh)  _ 
tan  ( ph) 

where 

CO2 

p2  =  —  —  h2 ,  c[2  — 
c2 

and  h,  k,  cl,  ct,  cp,  and  co  are  the  half-plate  thickness,  wave  number,  velocities  of  longitudinal 
and  transverse  modes,  phase  velocity,  and  wave  circular  frequency,  respectively. 

The  phase  velocity  is  found  numerically  by  solving  the  real  roots  of  the  dispersion  equation  as  a 
function  of  material  properties,  frequency,  and  material  thickness.  Group  velocity  dispersion 
curves  are  found  from  the  phase  velocity  (9). 

Phase  and  group  velocity  dispersion  curves  for  a  2.0-mm-thick  aluminum  plate  are  shown  in 
figure  1,  where  longitudinal  velocity  (cl),  transverse  velocity  (ct),  and  density  (p)  of  the  plate  are 

■j 

equal  to  6420  m/s,  3040  m/s,  and  2700  kg/m  ,  respectively.  These  curves  were  generated  and 
plotted  using  a  commercial  software  package  (10). 
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The  figure  shows  the  first  three  symmetric  and  antisymmetric  mode  dispersion  curves.  At  low 
frequencies,  the  wave  velocity  of  the  first  symmetric  mode  (So)  is  nearly  nondispersive,  while  the 
wave  velocity  of  the  first  antisymmetric  mode  (Ao)  is  highly  dispersive.  At  higher  frequencies, 
phase  velocity  of  both  zero-order  modes  approaches  the  Rayleigh  wave  velocity,  cr.  which  has  a 
value  of  ~2.9  km/s  ( 6)  for  aluminum  with  a  Poisson’s  ratio  v  =  0.33. 


(a) 


(b) 


Figure  1.  Dispersion  curves  for  2. 0-mm-thick  aluminum  plate 
(a)  phase  velocity  and  (b)  group  velocity  vs. 
frequency,  where  SO  and  AO  are  first-order  symmetric 
and  antisymmetric  modes,  respectively. 
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3.  Finite-Difference  Model  and  Simulations 


A  finite-difference  model  was  used  to  simulate  wave  propagation  in  a  thin  aluminum  plate  with 
notch.  For  this  modeling  and  simulation,  commercial  software  (11)  was  used  based  on  a 
published  algorithm  (12)  that  solves  a  2-D  (plane  strain)  acoustic  wave  equation.  A  schematic 
diagram  of  the  finite-difference  time  domain  model  is  shown  in  figure  2.  In  the  model,  the 
ultrasonic  source  and  receiver  were  each  10  mm  in  diameter  and  in  planar  contact  with  the 
sample  surface  positioned  100  mm  apart.  The  aluminum  plate  was  2.0  mm  thick  and  600  mm 
wide.  A  wide  width  was  used  to  prevent  interference  of  edge  reflections.  Aluminum  properties 
of  the  material  used  in  the  model  were  longitudinal  velocity  of  6420  m/s,  transverse  velocity  of 
3040  m/s,  and  density  of  2700  kg/m  .  The  simulated  damage  (the  notch)  had  a  width  of  6.0  mm 
and  depth  from  the  sample  back  surface  that  ranged  from  0.0  to  1.0  mm  in  0.1 -mm  steps.  Notch 
depth  translated  to  a  depth- to-thickness  ratio  of  0%  to  50%.  To  vary  the  location  of  the  notch, 
the  source  and  receiver  transducers  were  repositioned  but  always  kept  100  mm  apart.  For 
example,  the  notch  was  centered  at  0.0  along  the  axis,  the  source  was  positioned  at  -40  mm,  and 
the  receiver  positioned  at  60  mm. 
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Figure  2.  Finite-difference  model  schematic. 

Two  different  transmitter  or  source  signals  were  used.  A  sine  Gaussian  pulse  and  sine 
exponential  pulse  were  available  as  functions  in  the  software.  Source  signals  and  spectrums  for 
the  sine  Gaussian  pulse  and  sine  exponential  pulse  are  shown  in  figures  3  and  4,  respectively. 
The  sine  Gaussian  pulse  spectrum  has  a  full-width-half-height  (FWHH)  bandwidth  of  -0.32 
MHz  and  ranges  from  0.0  to  1 .0  MHz.  The  sine  exponential  pulse  spectrum  has  a  FWHH 
bandwidth  of -0.43  MHz  and  frequency  range  from  0.0  to  -2.0  MHz,  with  low  amplitude 
frequency  components  above  1  MHz. 
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Figure  3.  Sine  Gaussian  function  and  its  spectrum  used  in  the  simulation 
as  a  source  function. 
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Figure  4.  Sine  exponential  function  and  its  spectrum  used  in  the 
simulation  as  a  source  function. 
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Example  temporal  snapshots  of  guided  waves  propagation  in  the  sample  with  a  centrally  located 
notch  are  shown  in  figure  5.  This  figure  shows  how  the  wave  spreads  as  a  function  of  time.  An 
example  receiver  signal  is  shown  in  figure  6,  where  figure  6a  shows  a  signal  for  no  damage  (or 
baseline)  and  figure  6b  shows  a  signal  for  a  centrally  located  notch  with  a  0.5  mm  deep.  For 
these  signals,  the  So  mode  is  located  at  ~20  jus  and  the  Ao  mode  starts  at  ~30  ps.  A  series  of 
signals  were  generated  for  various  notch  depth  and  position  and  were  examined  in  time, 
frequency,  and  time-frequency  domains. 


Source  Receiver 


Guided  wave  mode:  A0  S0 


Figure  5.  Typical  snapshots  of  the  guided  wave  propagating  in  the  model. 


Figure  6.  Simulated  time-domain  signal  showing  the  primary  symmetric  and 
antisymmetric  modes  for  (a)  baseline  signal  without  damage  and 
(b)  a  notch  with  0.5-mm  depth.  (Source  signal  was  a  sine 
exponential  pulse.) 
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4.  Time  and  Frequency  Analysis  and  Results 


Time  and  frequency  domain  analysis  of  signals  similar  to  the  ones  shown  in  figure  6  was 
performed.  In  this  analysis,  wave  transmission  coefficients,  Fourier- transform  magnitude, 
spectral  energy,  and  average  powered  vs.  notch  depth  were  considered.  The  A0  and  So  mode 
transmission  coefficients  were  obtained  by  dividing  the  subject  mode  amplitude  by  its 
corresponding  baseline  amplitude.  The  Ao  mode  transmission  coefficient  over  a  range  of  notch 
depth  is  shown  in  figure  7a  and  shows  an  almost  constant  value  to  a  notch  depth  of  about  0.5  mm 
and  decreasing  values  for  deeper  notch  depths.  A  plot  of  transmission  coefficients  for  the  So 
mode  in  figure  7b  shows  a  similar  result.  Spectral  analysis  for  a  series  of  notch  depths  is  shown 
in  figure  8.  The  Fourier-transform  magnitude  is  shown  in  figure  8a.  Overall,  the  magnitude 
trend  is  difficult  to  determine,  but,  in  general,  there  is  a  decreasing  magnitude  for  increasing 
notch  depth.  The  spectral  energy  in  figure  8b  also  shows  decreasing  energy  for  increasing  notch 
depth.  Average  signal  power  as  a  function  of  notch  depth  is  shown  in  figure  9.  This  shows  a 
similar  trend  as  depicted  in  figures  7-9.  These  plots  showed  time  or  frequency  aspect  of  a 
guided  wave  signal  vs.  notch  depth.  They  showed  a  flat  response  or  no  change  in  signal 
amplitude,  spectral  magnitude,  and  average  power  up  to  a  damage  depth  of  ~0.5  mm  (25%  of 
material  thickness),  making  them  insufficient  to  determine  shallow  notch  depth. 
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Figure  7.  Transmission  coefficient  for  (a)  A0  mode  vs.  notch  depth  and  (b)  So 
mode  vs.  notch  depth. 
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Figure  8.  Spectral  analysis  of  signal  shows  (a)  the  Fourier-transform  magnitude 
vs.  notch  depth  and  (b)  spectral  energy  vs.  notch  depth. 
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Figure  9.  Average  power  vs.  notch  depth. 
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5.  Integrated  Time-Frequency  Analysis  and  Results 


A  time-frequency  representation  (TFR)  (73)  is  another  method  used  for  signal  analysis.  It  has 
been  used  to  identify  Lamb  wave  modes  and  dispersion  curves  which,  in  turn,  can  be  used  to 
calculate  material  properties.  The  TFR  maps  a  one-dimensional  time-domain  signal  into  a  2-D 
function  of  time  and  frequency  displayed  as  a  gray-scale  or  false-color  image.  The  TFR  used 
here  is  the  short-time  Fourier  transform  (STFT)  or  its  magnitude  squared  the  spectrogram  and  the 
reassigned  spectrogram  (14-19). 

For  reference,  the  STFT  is  as  follows  (20): 

/< oo 

x(t)  co(t  -t)  e-j°“  dt ,  (4) 

-OO 

where  co(f)  is  the  window  function,  typically  a  Hann,  Hamming,  or  Gaussian  window,  and  x(t)  is 
the  signal  to  be  transformed.  The  spectrogram,  the  STFT  magnitude  squared,  is  as  follows  (20): 

spectrogram  {x(t)}  =  | X(x,co)|" .  (5) 

Spectrogram  resolution  depends  on  the  window  size,  type,  and  overlap.  The  width  of  the 
window  relates  to  frequency  resolution  (the  separation  of  frequency  components)  or  time 
resolution  (the  time  at  which  frequencies  change).  A  wide  window  gives  better  frequency 
resolution  but  poor  time  resolution.  A  narrow  window  gives  good  time  resolution  but  poor 
frequency  resolution.  As  a  result,  a  tradeoff  exists  between  improved  time  and  frequency 
resolution.  The  reassigned  spectrogram  sharpens  time  and  frequency  estimates  and  thus  can  be 
considered  a  postprocessing  on  the  STFT  distribution  for  improve  readability.  One  reassignment 
method  (18,  19,  21)  uses  a  moving  window  to  maximize  the  coefficients  in  a  local  time- 
frequency  distribution,  and  another  method  (22)  uses  results  of  two  STFTs  with  slightly  modified 
window  functions. 

The  signals  in  figure  6  with  and  without  damage  (notch  depth  of  0.5  mm  or  25%  of  material 
thickness)  were  mapped  to  the  time-frequency  domain.  The  spectrogram  for  the  baseline  signal 
in  figure  6a  is  shown  in  figure  10a.  In  this  spectrogram,  the  dispersion  group  velocity  curves  the 
So,  Si,  A0,  and  Ai  modes)  from  figure  lb  are  superimposed  over  the  time-frequency  map  and 
illustrate  the  time-frequency  and  group  velocity  mapping  for  the  given  source  bandwidth. 
Dispersion  curve  group  velocity  was  converted  to  time  using  the  source  receiver  separation 
distance  as  a  conversion  parameter.  The  intensity  map  shows  the  primary  Ao  and  So  Lamb  wave 
modes  and  some  higher  modes,  with  the  So  mode  arriving  at  about  18  ps  and  the  A0  mode 
arriving  at  about  30  ps.  Frequency  dispersion  of  the  Ao  can  be  seen  with  the  higher  frequency 
modes  arriving  before  the  lower  frequency  modes.  The  spectrogram  for  the  signal  with  damage 
is  shown  in  figure  10b.  This  figure  is  similar  to  figure  10a  except  that  the  region  between  the  So 
and  Ao  modes  at  about  25  ps  shows  a  mode  conversion  resulting  from  the  notch. 
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Figure  10.  Spectrogram  for  (a)  baseline  signal  and  (b)  signal  with  0.5-mm-deep  notch  for  a  sine 
exponential  source  signal. 

To  see  the  result  of  a  different  source  signal,  simulated  signals  were  generated  for  the  same 
geometry  as  the  ones  shown  in  figure  6,  but  with  a  sine  Gaussian  source  function.  The  result 
mapped  to  the  time-frequency  domain  in  figure  1 1  shows  a  reduced  bandwidth  as  a  result  of  the 
source  signal.  As  in  figure  10,  figure  1 1  shows  the  mode  conversion  from  the  notch  in  the  area 
between  the  So  and  A0  modes. 

When  a  reassigned  spectrogram  was  applied,  the  TFR  coefficients  in  figures  10  and  1 1  follow 
the  dispersion  curves  along  a  tighter  focus.  These  results  are  shown  in  figures  12  and  13  for  the 
sine  exponential  and  sine  Gaussian  pulse,  respectively.  Although  these  results  show  a  tighter 
focus,  the  reassignment  process  dissipates  and  the  mode  converts  the  signal,  making  the 
reassignment  spectrogram  less  useful  for  the  current  purpose. 


Time  (psec) 

(a) 


Time  (psec) 

(b) 


Figure  11.  Spectrogram  for  (a)  baseline  signal  and  (b)  signal  with  0.5-mm-deep  notch  for  a  sine  Gaussian 
source  signal. 
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Figure  12.  Reassigned  spectrogram  for  (a)  baseline  signal  and  (b)  signal  with  0.5-mm-deep  notch  for  a  sine 
exponential  source  signal. 


Figure  13.  Reassigned  spectrogram  for  (a)  baseline  signal  and  (b)  signal  with  0.5-mm-deep  notch  for  a  sine 
Gaussian  source  signal. 

Damage  or  notch  sensitivity  was  examined  through  the  spectrogram  by  processing  guided  wave 
signals  generated  for  a  series  of  notch  depths  and  the  sine  exponential  source  signal,  which  has  a 
larger  bandwidth  than  the  sine  Gaussian  source  signal.  Of  the  smaller  notched  depths  examined 
0.05  and  0.10  mm  deep,  only  the  latter  showed  the  effect  of  the  mode  conversion.  The 
spectrogram  for  the  0. 10-mm  notch  is  shown  in  figure  14a,  and  a  0.9-mm  notch  is  shown  in 
figure  14b.  The  result  is  similar  to  figure  10b,  with  the  mode  converted  signal  appearing  at  about 
25  ps.  Comparing  the  mode-converted  signal  in  figures  10b  and  14,  it  appears  that  the  intensity 
of  the  mode-converted  signal  is  increasing.  However,  when  the  intensity  was  examined  for  the 
series  of  notch  depths,  the  mode-converted  signal  intensity  did  not  consistently  increase  as  depth 
increased.  This  was  most  likely  due  to  wave  addition  and  subtraction  and,  thus,  varying  wave 
intensity. 
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Figure  14.  Spectrogram  for  a  signal  with  (a)  0.1-mm  and  (b)  0.9-mm-deep  notch  and  sine  exponential 
source  signal. 

The  damage  or  notch  for  the  spectrogram  images  was  centered  between  the  source  and  receiver. 
When  the  notch  was  offset  from  the  center,  the  temporal  location  of  the  mode  converted  signal 
shifted.  Figure  15  shows  spectrograms  for  a  0.1-mm-deep  notch  centered  in  figure  15a,  offset 
10  mm  in  figure  15b,  and  offset  20  mm  in  figure  15c.  In  this  sequence,  the  temporal  location  of 
the  mode  converted  signal  in  figure  15a-c  moves  from  ~25,  to  27,  to  29  jas,  respectively.  These 
temporal  locations  are  difficult  to  determine  accurately  because  the  spectrogram’s  So  and  Ao 
mode  components  are  broad  and  fuzzy. 


Figure  15.  Damage  indicator  temporal  location  as  a  function  of  a  0.1-mm-deep  notch  physical  location 
between  the  source  and  receiver  (a)  centered  or  no  offset,  (b)  10-mm  offset,  and  (c)  20-mm 
offset,  with  damage  temporal  locations  of  ~25,  27,  and  29  us,  respectively. 
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6.  Conclusions 


Time,  frequency,  and  time-frequency  analysis  methods  of  ultrasonic  guided  waves  for  detecting 
a  small  notch  were  present.  Ultrasonic-guided  wave  signals  were  numerically  generated  using  a 
finite-difference  model.  In  this  model,  damage  was  simulated  in  a  rectangular  notch  of  varying 
depth  and  position  between  a  source  and  receiver  transducer.  The  signals  displayed  the  classical 
symmetric  and  antisymmetric  modes.  Analysis  of  these  signals  as  a  function  of  notch  depth 
showed  a  flat  response  or  no  change  in  signal  amplitude,  spectral  magnitude,  and  average  power 
up  to  a  damage  depth  of  ~0.5  mm  (or  25%  of  the  sample  thickness).  For  greater  damage  or 
notch  depths,  larger  changes  were  observed.  This  result  showed  that  gross  damage  could  be 
detected,  but  not  shallow  damage  less  than  25%  of  material  thickness. 

Time-frequency  representation  of  signals  was  more  successful  in  displaying  signal  changes  due 
to  shallow  depth  damage.  Spectrogram  images  displayed  artifacts  from  the  notch  for  shallow 
depths  as  small  as  0.1  mm  (or  5%  of  sample  thickness).  Of  the  two  source  waveforms  used,  the 
sine  exponential  pulse  showed  a  better  localization  of  signal  artifacts  from  the  notch,  although 
this  may  be  an  effect  of  image  intensity  settings.  The  reassigned  spectrogram  focused  the  time- 
frequency  data  tighter  along  the  dispersion  curves  but  seemed  to  dissipate  the  data  from  the 
notch  as  a  result  of  the  reassignment  algorithm;  this  was  not  beneficial  to  the  analysis. 
Spectrogram  sensitivity  to  the  notch  position  was  examined.  This  illustrated  that  the  notch  could 
be  tracked  to  a  temporal  position  between  the  source  and  receiver,  but  accuracy  was  poor  due  to 
the  cloud  nature  of  the  TFR. 

Overall,  the  time-frequency  representation  displayed  signal  variations  due  to  a  shallow  notch 
depth  as  small  as  0. 1  mm,  5%  of  material  thickness  which  is  5X  smaller  than  what  could  be 
detected  by  time  and  frequency  analysis  methods  alone.  Thus,  TFR  methods  can  be  used  to 
interrogate  a  structure  and  detect  small  damage  which  can  then  be  repaired  to  prevent  serious 
damage  to  equipment  and  personnel. 
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